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Abstract 

Boost-invariant (1+1) dimensional solutions for the Disoriented Chiral Con- 
densate (DCC) are obtained numerically, in the context of the SU(2)l (g) 
SU(2)r non-linear sigma model at 0{p i ) in the momentum expansion. We 
ignore the mass terms in the Lagrangian, as we are mainly interested in the 
behavior of the solutions for small values of proper time r. The solutions ob- 
tained at 0(p 4 ) are matched to those of 0(p 2 ) at a late proper time r 3> 1/77%, 
where m n = 140 MeV is the mass of the pion. We find that at C(p 4 ) the 
solutions for the DCC do not have singular behavior at early proper times 
r <C l/m n . The solutions indicate that for r < (0.5 — 0.8) fm the 0(p 4 ) 
corrections become important. We take the sizes of the field derivatives to 
be indicators of the validity of the momentum expansion. Thus, we deduce 
that the C(p 4 ) solutions can be used to represent the qualitative behavior of 
the DCC down to proper times of about 0.2 fm. Since below ~ 0.2 fm the 
formalism is not reliable, we conclude that the inclusion of higher order terms 
beyond 0(p 4 ) is not needed to extend the validity of the solutions to earlier 
proper times. 
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I. INTRODUCTION 



In quantum chromodynamics (QCD), the Lagrangian for strong interactions is invariant 
under global SU (Nf) L ®SU (Nf) R chiral transformations, where Nf is the number of massless 
flavors. However, the non-zero masses of quarks in the Standard Model explicitly break this 
global symmetry. In addition, chiral symmetry is dynamically broken at a scale A x ~ 1 GeV 
where the vacuum expectation value of the quark bilinear (q % qj) becomes non-zero. Thus, 
for Nf massless quarks and at energies above A x , the SU(Nf)i ® SU(Nf)n symmetry of 
strong interactions is restored. The restoration of chiral symmetry at non-zero temperatures 
has been studied in the context of the theory of phase transitions and critical phenomena 
IJ. According to universality arguments and numerical simulations, in order for the chiral 
transition to be a second order phase transition, only two quark flavors, namely the up and 
the down quarks, can be treated as massless 0. Since the up and the down quark masses 
are nearly zero compared to the scale of symmetry breaking, the global SU(2)i <g> SU(2)r 
is an approximate symmetry of QCD at energy scales above A x . 

If the masses of the quarks were zero the condensate (q i qj) would not have any preferred 
direction in the vacuum, under SU{2)l ® SU(2)r. However, the small but non-zero quark 
masses provide the QCD vacuum with a direction. It has been argued that in certain high 
energy collisions, such as high energy pp |2j or relativistic heavy ion collisions domains 
of chiral condensates that do not point in the direction of the QCD vacuum may form and 
grow to a volume of (a few fermi) 3 . These domains are referred to as Disoriented Chiral 
Condensates (DCC's). The eventual decay of DCC's into a large number of pions is predicted 
to have a distinct experimental signature, namely large fluctuations in the ratio / of neutral 
pions to the total number of pions. 

The non-linear sigma model, based on SU{2)l <g> SU(2)r, has the essential features 
necessary for describing low energy QCD phenomena, and has been used to describe the 
evolution of the DCC after its formation [0-^,0-0. In Refs. |5|f|, the authors assume that 
a relativistic collision can be described by two thin infinite slabs, representing the highly 
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Lorentz-contracted hadrons or nuclei, that collide at the center of mass, and continue along 
the beam axis in opposite directions. The DCC forms in the region of spacetime between 
the receding hadrons or nuclei. These authors find boost-invariant (1+1) dimensional clas- 
sical pion field solutions of the non-linear sigma model at the leading order in a derivative 
expansion, without the mass term. The solutions they obtain exhibit violent oscillations for 
small values of the proper time r, which can be taken as a sign of the breakdown of the 
formalism at early proper times, where higher order effects become important. 

In this paper, we include the four-derivative terms, ignoring the mass terms, in the 
Lagrangian for the £77(2) £, (g> SU(2)r non-linear sigma model. We numerically solve the 
Euler-Lagrange equations to obtain the DCC solutions which are shown to be non-oscillatory 
at early proper times, r <C l/m^, and we present an analytic explanation for this behavior. 
We also show that the corrections from terms with four derivatives are important at proper 
times earlier than about (0.5-0.8) fm. The 0(p 4 ) solutions for the fields parametrizing the 
DCC have derivatives that diverge for r — > 0, as in the case of the previously obtained 0(p 2 ) 
solutions. However, the proper times at which the magnitudes of the derivatives are deemed 
too large for a reliable momentum expansion are smaller than those of the leading order 
solutions. In the rest of this paper, by 0(p 2 ) solutions we mean those obtained from the 
Lagrangian with terms that have at most two derivatives. The solutions that are obtained 
from the Lagrangian that contains terms quadratic and quartic in derivatives are referred 
to as the 0(p 4 ) solutions. 

In the next section, we outline the formalism of the non-linear sigma model at 0(p A ). 
In Section III, we present the 0{p 2 ) solutions obtained in Refs. |5]||. Section IV, where we 
present our numerical solutions, includes a discussion of our choice of boundary conditions 
and a comparison of our solutions with those of the leading order. Section V contains a 
summary of our results and some concluding remarks. 
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II. THE NON-LINEAR SIGMA MODEL AT 0(p A ) 



In this section, we establish the formalism used in this paper to study the evolution of 
the DCC at the next to leading order in momentum expansion. We use notation similar to 
that of Ref . [0] . To represent the effects of chiral symmetry breaking and the pion fields, we 
introduce the fields cr(x) and n(x), respectively. We represent the pion field by 

tt(x) = fn(x) sin 9(x), (2.1) 

where at leading order / is the pion decay constant f n = 93 MeV, n is a unit isovector field 
\n\ 2 = 1, and the field 9 is an angle. The isovector n determines the orientation of the pion 
field in isospace. The a and 7? fields are related by a 2 + |7?| 2 = f 2 . We define the field £ by 

f n E = a + tf-rr, (2.2) 

where r % for i = 1, 2, 3 are the Pauli matrices. 

The Lagrangian for the non-linear sigma model at 0(p 2 ), without the mass terms, can 
then be written as 

f 2 

£(2) = i^Tr(d^d^). (2.3) 

In terms of the 9 and the n fields, we have 

£(2) = 11^98^9 + sin 2 6 d^n ■ d»n) + ^-(\n\ 2 - 1), (2.4) 

where A is a Lagrange multiplier. 

To go to a higher order in the momentum expansion, 0(p A ), we should include opera- 
tors with four derivatives, or two derivatives and one insertion of the quark mass matrix. 
However, we continue to ignore the mass terms, since we are mainly interested in the early 
evolution of the DCC, which corresponds to regions of large momenta. The 0(p A ) contribu- 
tion to the Lagrangian of the system is then given by 

£( 4 ) = ft \Tr(d lt Ed»E i )] 2 + /3 2 Tr(a M E^E t )Tr(^E^E t ), (2.5) 



where, in the notation of Ref. |10| , fa = a\ + fa = «2, and aj, i = 1,2,3, are the 

coefficients of the 0(p 4 ) terms for the chiral Lagrangian under SU(3)l ® SU(3)r |10[. In 
terms of 9 and n, we have for 



£^=4fa \(d, 



+ 2 sin 2 9 {8^9) (d u n ■ d v n) + sin 4 9 {d^n ■ d^n) 



+4 fa {d^d^f + sin 4 9 (<9 M n • d u n) {d"n ■ d y n) + 2 sin 2 9 djdj (d»n ■ d v n) 



(2.6) 



The coefficients fa$ get renormalized by one-loop diagrams coming from the 0(p 2 ) La- 
grangian, including the mass terms, and are phenomenologically determined. In this work, 
we are only interested in the classical behavior of the pion field, and hence we do not consider 
the loop effects. As our results represent the qualitative behavior of the DCC, we only give 
an order of magnitude estimate for the typical size \/3\ of fa and fa, and ignore the quantum 
corrections coming from loop diagrams. We expect this approximation to be qualitatively 
valid, since by large N c arguments, where N c is the number of colors, the effects of loop 
corrections are suppressed by iV" 1 . 

In order to estimate the fa^, we demand, in a systematic momentum expansion of the 
Lagrangian 



£ P = (^-)p 2 + \P\p 4 + o(p 6 ) 



(2.7) 



that all terms be of the same order at the chiral symmetry breaking scale A x . That is, 



and hence 
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f 2 

\ ~4 " \ ' 



f 2 

J 7T 



4A 2 ' 



(2.8) 



For f w = 93 MeV and A x = 1 GeV, we get 



2 x 10" 



(2.9) 



For phenomenologically relevant energy scales, fa + fa > [ID|. Since in this paper it is the 
combination fa + fa that appears in the equations of motion, we take (3 = fa + fa > 0. 



III. THE 0(p 2 ) SOLUTIONS 



The 0(p 2 ) Euler-Lagrange equations of motion derived from Eq. fl2.4j ) are given by 

DO = sm9 cos Od^n- d^n (3.1) 

and 

d^m 2 9d^n) = \n. (3.2) 

In Refs. [0,0, the authors have assumed that the DCC solutions for a relativistic hadronic 
or nuclear collision have transverse symmetry with respect to the beam direction, where the 
colliding particles are idealized as two highly Lorentz-contracted slabs of infinite transverse 
extent. This idealization makes the problem (1+1) dimensional. We take these dimensions 
to be time t and the beam direction x. With the further condition of boost-invariance, the 
DCC solutions [||,|| become functions of only proper time r = \/t 2 — x 2 , that is n = n(r), 
and 9 = 9{r). Note that for a function = <f>{r) we have <9^0 = (x M /r)0', where a prime 
denotes a derivative with respect to r; 0' = dip /dr. Here, we simply mention the solutions 
to Eqs.( p.l|) and ( |3.2|) ; the details of the solution are found in Refs. [|]||. 
The angle 9 is given by 

cos#(r) = (&/«)cos[k1ii(t/to) +0 O ], (3.3) 



where k = ^/ a 2 + 6 2 , a = \a\, and b = \b\; a and b are two arbitrary constant vectors. In 
addition, cos^o — i K /b) cos^(ro), and tq is some arbitrary proper time. The choice of the 
coordinates in isospace is such that a, c = b x a, and b define a right-handed coordinate 
system, and n a = 0. We also have 



n b (r) = (~). 



K 2 



(3.4) 



b J \a 2 sin 2 9(t)' 

The solution for n c (r) is obtained from the constraint \n\ 2 = 1. In this work, we take 
cos#(to) = implying sin 6^ = [§, which will make the above solutions (|3.3|) and ( p.4|) 
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of the same form as those presented in Ref. ||. Note that for a = 0, we have n = Constant 
as the solution. In the rest of this work, we refer to n a , n b , and n c by nj, n 2 , and n 3 , 
respectively. 

The solution for 9(t) given by Eq.( |3.3| ), oscillates rapidly near r = 0. Small values of r 
correspond to regions in spacetime that are close to the highly energetic nuclei. We expect 
the theory to break down for small r, corresponding to large momenta, since our theory is 
valid only at low momenta. Thus, it may be tempting to interpret tq as a typical proper 
time below which the theory becomes unreliable ||, as signaled by the "rapid" oscillations. 
However, the onset of these "rapid" oscillations is scale-dependent, and cannot reliably 
determine r , in the above sense. 



IV. THE CO 4 ) SOLUTIONS 

In this section, we will show that the solutions to the 0(p 4 ) equations show no divergent 
behavior near r = 0. This does not mean that we can trust the qualitative behavior of 
the solutions for arbitrarily small r. Instead, we note that our theory is an expansion 
in derivatives, and thus we will take the magnitudes of the derivatives of the pion field 
parameters 9 and n to be better indicators of the range of the validity of our solutions. We 
will thus demonstrate that the 0(p 4 ) solutions stay reliable down to a length scale of about 
0.2 fm. 

The Euler-Lagrange equations of motion for 9 and n are derived from the Lagrangian £, 
given by 

£ = £ (2) +£ (4) , (4.1) 
where £ includes the derivative terms up to 0(p 4 ). The equation for 9 is 



fl ^ + fl"-sin#cos#K| 2 J +16(/3i + /32)[(~ + 0"J (9' 2 + sin 2 9\n'\ 2 ) 
+29' 2 9" + 2 sin 2 9 9'{n ■ n") + sin 9 cos 9\n\ 2 (9' 2 - sin 2 9\n'\ 



'9' 



(4.2) 
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For n, we get the following equation 



J 7T 



2 cos 9 9' fv + sin 9 I h n " + \n'\ z n 



+ 16(/3 1 + /3 2 K sin9 



29'9"n' 



-9' 2 (- + n") + 4 sin 9 cos 9 9' '\n' f n' + sin 2 9\n'\ 2 f-+rf 



+2sitf0(n'-n")n' + (sitf 9\n\ z + 9' ) \n'\ z n 



T 



+ 2cos66'n'} = 0, (4.3) 



where we have used the constraint n ■ n" = — \n'\ 2 . 

To solve the above coupled non-linear ordinary differential equations (}4.2|) and (|4.3| ), we 
need to specify the boundary conditions. In this work, we choose the boundary conditions 
for the 0(p 4 ) fields 9 and n and their derivatives to be the values of the corresponding 
0{p 2 ) solutions evaluated at a "late" proper time 77, where the 0{j> 2 ) and 0{p A ) solutions 
approximately coincide. We choose To = \j2m^^ as a typical proper time where we expect 
the higher order interactions to become important, j3i + = (3 = 2 x 10~ 3 , and a = b = 1. 
We pick T\ = IOtq. In this way, the solutions of the 0(p 2 ) equations have definite values at 
77. Numerically, we have r m 0.7 fm, and 77 ~ 7 fm, at which we expect the 0{p 2 ) solutions 
to approximate the 0(p A ) solutions with good accuracy. 

In Figs.(l) through (6), we present our numerical solutions for the fields 9 and n and 
their derivatives. Note that as a result of current conservation relations [§], n\ — and 
n[ = 0, as in the case of the leading order solutions. Each figure contains the 0(p 2 ) solution, 
represented by the dashed line, and the 0(p 4 ) solution, represented by the solid line. The 
0{p 2 ) solutions are obtained numerically and agree with those of the analytic expressions 
in Eqs.Q and Q . 

The solutions presented here are computed for r 6 [lO" 4 ^/???,,,-] MeV -1 . Figures (1), 
(2), and (3) show that the rapid oscillations of the 0(p 2 ) solutions for 9 and ft, near r = 0, 
no longer arise in the 0(p A ) solutions, where the inclusion of the higher order terms seems 
to stabilize the solutions. We tested the stabilization of the solutions by comparing the 
values of the fields at proper times r < 10~ 4 MeV -1 , and observing that whereas the leading 
order solutions oscillate more rapidly for smaller r, the values of the 0(p 4 ) solutions do not 



oscillate or change significantly. 

The 0(p 4 ) and the 0(p 2 ) field derivatives presented in Figs. (4), (5), and (6) show di- 
vergent behavior near r = 0. The analytic 0(p 2 ) solutions of Eqs.( |3.3| ) and ( |3.4| ) yield 
divergent and oscillatory derivatives for r = 0. To check the divergence of the 0(p A ) deriva- 
tives, we examined their behavior for r < 10~ 4 MeV -1 . We found that the magnitudes of 
the derivatives did not stabilize and continued to grow with decreasing r. 

In order to understand the behavior of the 0(p 4 ) solutions mentioned above, we examine 
Eq.( [4.2|) in the limit r — > 0, for the case of constant n. The solutions with constant n are 
related to those with spacetime dependence by chiral rotations ||. In this case, Eq.( |1.2| ) 
reduces to 

fl + + 16(fa + fa) + 39"^ 9' 2 = 0. (4.4) 
Let us assume that for small r, the behavior of the field 9 is given by 

o=~e+{l)\ (4.5) 

where 9 and f are constants. We expect to find a solution for 9 with < p < 1 that tends 
to a constant 9 but has a divergent derivative, as r — > 0. Upon substituting the expression 
in Eq. ( |4.5[) for 9 into Eq. fl4.4j) , we get the following equation 

fl[p + p(p ~ 1)] + 16(A + fo)\p + MP ~ W (Vf P Q * = 0. (4.6) 

In the absence of the 0{p i ) terms, we only have the terms proportional to f%. In this case, we 
get p = 0, which is consistent with the solution 9 ~ ln(r/f) of the equation 6'/t + 6" = 0. 
However, we see that for small r, the terms proportional to fl\ + fa will be dominant if 
< p < 1. Thus, for t —>■ 0, we ignore the terms proportional to f%. We then get 

P = | (4-7) 

in agreement with our prior assumption that < p < 1. 

The numerical solutions presented here are obtained in the context of a systematic deriva- 
tive expansion. Thus, it is the magnitudes of the derivatives that establish the region of 
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validity of the solutions. We take the maximum magnitude of the derivative of a field below 
which the expansion is valid to be p max ~ 500 MeV. The graphs in Figs. (4), (5), and (6) 
show that the magnitudes of the 0{p A ) field derivatives stay below p max for values of r 
down to t = 10" 3 MeV" 1 ~ A" 1 , which is as small a proper time as we can consider in 
our formalism. In contrast, the magnitudes of the 0(p 2 ) field derivatives begin to exceed 
Pmax at proper times r < (2 — 4) x 10~ 3 MeV -1 , indicating a smaller range of validity for 
these solutions. Hence, our results suggest that the 0(p i ) solutions can be used to study 
the qualitative evolution of the DCC down to a length scale of ~ 0.2 fm, but the 0(p 2 ) 
solutions lose their qualitative validity below a length scale of ~ (0.5 — 0.8) fm. 



V. CONCLUSIONS 

We derived the 0(p 4 ) equations of motion for the DCC produced in an idealized rel- 
ativists collision, using the non-linear sigma model Lagrangian without the mass terms. 
We presented our numerical solutions for the C(p 4 ) equations of motion. These higher or- 
der corrections to the 0{p 2 ) solutions are only important for small values of proper time 
t <C l/m n . Hence, the absence of the mass terms, important only for r > l/m n , does not 
introduce significant qualitative changes in the early proper time solutions. Our 0(p 4 ) so- 
lutions for the fields 9 and n, parametrizing the DCC field configuration, stabilize for small 
values of proper time, whereas the 0{p 2 ) solutions oscillate rapidly with decreasing r. We 
presented a qualitative explanation for the behavior of the solutions by studying the equa- 
tions of motion, in a special case and in the limit where r — > 0. A measure of the validity of 
our derivative expansion is the magnitude of the field derivatives. We took p max ~ 500 MeV 
to be the maximum value for the magnitude of a field derivative beyond which the solutions 
cannot be trusted. Using this criterion, the 0(p 4 ) solutions were deemed reliable down to 
r ~ 0.2 fm, whereas the 0{p 2 ) solutions were considered no longer qualitatively valid below 
t ~ (0.5 — 0.8) fm. We take our solutions to represent a qualitative measure of the behav- 
ior of the DCC. Our results suggest that this qualitative behavior can be studied within 
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the non-linear sigma model down to a length scale of ~ 0.2 fm, once the 0(p 4 ) derivative 
corrections to the Lagrangian are included. A length scale of ~ 0.2 fm corresponds to the 
energy scale A x ~ 1 GeV, where chiral symmetry is restored and the non-linear sigma model 
formalism is no longer valid. Hence, we do not expect the inclusion of the higher order terms 
beyond 0{p A ) in the chiral Lagrangian to enable us to study even earlier proper times in 
the evolution of the DCC. 
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FIGURES 




0.005 0.01 0.015 0.02 0.025 0.03 0.035 
T [MeV -1 ] 

FIG. 1. The dashed and the solid lines represent the 0(p 2 ) and the C(p 4 ) solutions for 9, 
respectively. The 0(p 2 ) solution for 9 oscillates increasingly rapidly with decreasing r, as Eq.( |3.3| ) 
implies, while the 0(p 4 ) solution does not oscillate as r — > 0. Similar comments apply to the 
solutions for rti and 713 in Figs. (2) and (3), respectively. 




0.005 0.01 0.015 0.02 0.025 0.03 0.035 
T [MeV -1 ] 

FIG. 2. The dashed and the solid lines represent the 0(p 2 ) and the 0(p 4 ) solutions for rt2, 
respectively. 
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0.005 0.01 0.015 0.02 0.025 0.03 0.035 
X [MeV -1 ] 

FIG. 3. The dashed and the solid lines represent the 0(p 2 ) and the 0(p 4 ) solutions for ra 3 , 
respectively. 




0.002 0.004 



0.006 
x [MeV -1 ] 



0.008 



0.01 



FIG. 4. The dashed and the solid lines represent the 0(p 2 ) and the 0(p 4 ) solutions for 9', 
respectively. The 0(p 2 ) solution for 9' oscillates and has a divergent magnitude, while the 0(p 4 ) 
solution for 9' diverges in magnitude, but does not oscillate. Similar comments apply to the 
solutions for n' 2 and n' 3 in Figs. (5) and (6), respectively. 
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FIG. 5. The dashed and the solid lines represent the 0(p 2 ) and the 0(p 4 ) solutions for n' 2 , 
respectively. 



1000 




0.002 0.004 0.006 0.008 0.01 
x [MeV -1 ] 

FIG. 6. The dashed and the solid lines represent the 0{p 2 ) and the 0(p 4 ) solutions for n' 3 , 
respectively. 
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